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Objectives 

The objective of this research is to develop analysis procedures to investigate the coupling of 
composite and smart materials to improve aeroelastic and vibratory response of aerospace 
structures. The structural modeling must account for arbitrarily thick geometries, embedded and 
surface bonded sensors and actuators and imperfections, such as delamination. Changes in the 
dynamic response due to the presence of smart materials and delaminations is investigated. 
Experiments are to be performed to validate the proposed mathematical model. 

Introduction 

Laminated composite structural elements with surface bonded/embedded sensors and actuators 
offer great potential for static and dynamic control. Essential to designing these advanced 
structures are accurate and efficient mathematical modeling techniques. Imperfections, such as 
delaminations, are common to most composites and must be included in the analysis procedure. 
Several mathematical models have appeared in the literature for the analysis of adaptive beams and 
plates. These include investigations based on the classical theory [1,2], which is limited to the 
analysis of thin plates, first order Mindlin type analyses [3,4] and potentially expensive layer-wise 
theories [5,6], A hybrid theory has also been presented by Mitchel and Reddy [7], Higher order 
refined theory is both accurate and computationally efficient and have been shown to be useful for 
modeling smart composite laminates [8,9], Finite element models [3,5,9] are attractive since they 
can include practical geometries and boundary conditions. 
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Several efforts for modeling delamination of composite laminates have been presented in the 
literature. Although three dimensional approaches [10,1 1] are more accurate than two dimensional 
theories [12-15], their implementation can be very expensive for practical applications. The layer- 
wise approach [16] is an alternative since it is capable of modeling displacement discontinuities. 
However, the computational effort increases with the number of plies. Recently, a refined higher 
order theory, developed by Gu and Chattopadhyay [17], was shown to be an efficient approach for 
modeling delaminations. This theory has been shown to agree well with both elasticity solutions 
[18] and experimental results [17]. 

Relatively little attention in the literature has been paid to modeling delaminated smart 
composites with surface bonded/embedded piezoelectric actuators and sensors [19,20]. Therefore, 
the goal of the current work is to investigate this important topic. The refined higher order theory 
is used. This theory allows for transverse shear deformation which is especially important in the 
analysis of composites thick composites. It also allows for both slipping and separation of the 
delaminated regions. The finite element implementation allows practical geometries and boundary 
conditions to be modeled as well. 


Analysis 

In the developed theory, the composite laminate is divided into several regions which represent 
the nondelaminated zone and portions of the composite laminate above and below the delamination. 
The refined theory is implemented in each region by satisfying the stress free boundary conditions 
which occur at the top and bottom surfaces of the laminate as well as at the delamination interface. 
This allows several of the higher order functions to be determined in terms of the lower order 
functions. Continuity conditions are formulated between the regions. These conditions are 
enforced in the finite element implementation using a penalty approach. The theory can be easily 
generalized to include multiple piezoelectric sensors, actuators and delaminations as shown in Fig. 
1 . 

The general displacement field with parameters indicated in Fig. 2 is defined as follows. 
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( 1 ) 


U(x,y,z) = u(x,y) + (z - c)[-^-w(x,y) + <J> x (x,y) 

V ox 

V(x,y, z) = v(x,y) + (z - c/~ w(x.y) + <t> y (x,y)^ 


+ (z - c) 2 u 2 (x, y) + (z - c) 3 u 3 (x, y) 
+ (z - c ) 2 v 2 (x, y) + (z - c ) 3 v,(x. y) 


W(x,y, z) = w(x, y) 


where U, V and W are the total displacements, u, v and w denote the midplane displacements of a 
point (x,y), the partial derivatives of w represent the rotations of normals to the midplane 
corresponding to the slope of the laminate and <j) x and <J> y represent the additional rotations due to 
shear deformation about the y and x axes, respectively. The quantities U2, U3, V2 and V3 represent 
higher order functions. This displacement field has the advantage of easily reducing to the well 
known classical theory if the higher order terms are eliminated. The thickness coordinate, z, is 
measured from the global midplane of the laminate and c is the local midplane where c = 0 for an 
undelaminated composite. 


Constitutive relations 

Equations which relate stress, strain, charge and electric field are derived from the electric 
enthalpy density function given as follows. 

^( £ ij’^‘) = "2 C ijkl £ ij £ kl — ~ kjjH j Hj (2) 

where Cjj andEj are components of the strain tensor and electric field vector, respectively andc,^, 
e ijk > ‘tnd ky are the elastic, piezoelectric and dielectric permittivity constants, respectively. The 
stress and charge are determined as follows. 

D - 8H 

D ' -air O) 

_ 

°' J = 3e~ ( 4 ) 

For an orthotropic composite laminate with piezoelectric layers that have orthorhombic mm2 
symmetry in the context of laminate theory, the constitutive relationships are simplified as follows. 
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where A are the induced strains (A] = A 2 = d 31 E 3 ). It is important to note that engineering normal 
(£1-3) ar *d shear (£4.6) strains are now used in the above equations, £3 is set to zero and only the 
inplane piezoelectric constants (d 3 j and d 32 ) are retained in the context of the current work. 


Refined displacement field 

To account for the effects due to delamination, it is necessary to partition the laminate into 
several different regions shown in Fig. 3 . These regions include the nondelaminated region (H), 
the region above the delamination (ft d2 ) and the region below the delamination (Q d3 ). The 

interface between the nondelaminated region and the delaminated regions, indicated by the dashed 
line in Fig. 3 , is denoted Qj. The general form of the higher order displacement field (Eqn. 1) is 
independently applied to each of these regions to describe displacements which account for 
slipping and separation due to the delamination. However, this displacement field does not 
necessarily satisfy the condition that the transverse shear stresses, 04 and G5, vanish at the top and 
bottom surfaces of the plate (z = ±h/2) as well as on the debonded interface surfaces (z = hj)in the 
delaminated region. That is, 

a 4 (x,y,±h / 2) = 0, a 5 (x,y,±h / 2) = 0 (7) 

o A {x,y,h l ) = 0 , a 5 (x,y,h,) = 0 (x,y)efl d (8) 

in which the superscripts “+” and denote the quantities related to the layers above and below the 
delamination, respectively, and refers to the region that contains the delamination. For 
orthotropic plates, these conditions are equivalent to the requirement that the corresponding strains 


4 



be zero on these surfaces. A refined displacement field is obtained by applying these boundary 
conditions for each region as follows. 


3w ; 


3x 

3W; 




3d: 


u 2 . - Vj +(z-Cj)J ^ L + < t , y i ]“ ( z-c i ) 3 ^~T$ yj (i = 1,2,3) 


3d, 


u 3 = Wj 


(9) 


where c ; = 


( a i +b i) 


is the local midplane and dj = bj - a, is the thickness of the region. The 


quantities aj and bj correspond to the bottom and top coordinates of the local region, respectively, 
as shown in Fig. 4. It is important to note that several of the higher order terms in the generalized 
displacement field are either found to be zero or defined in terms of lower order functions. The 

displacements in each region are defined by independent functions with parameters given as 
follows. 


when i = 1 


when i = 2 


when i = 3 



(x.y)eQ 

(x,y) s Q d2 
(x,y) e 


(10a) 

(10b) 

(10c) 


It is important to note that the thickness of the plate, h, may vary due to the presence of surface 
bonded sensors/actuators. 


Continuity Conditions 

Additional boundary conditions must be imposed to ensure continuity of displacements at the 
interface of the nondelaminated and delaminated regions (O,) as shown in Fig. 5. A vector of the 
displacements is constructed to simplify formulation of the boundary conditions as follows. 

'Uj ' 

U: = v, 

|wj Cl) 

The continuity conditions at the interface of the nondelaminated and delaminated regions («,) are 
imposed as follows. 
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u i=i =U 1=2 a 2 < z < min(b|,b 2 ) (x.y)en, (] 2 ) 

U i=1 =U i=3 max(a,,a 2 )<z<b, (x,y)en, ( 13 ) 

The above equations can be satisfied exactly for the classical theory since the displacement 
distribution through the thickness is linear. However, the displacement distribution using the 
refined theory is nonlinear and must be satisfied in an average sense as follows. An error function 
vector for the first of the above equations is formulated as follows. 

e = u i=i ~Uj=2 hi < z < h / 2 (x,y)er 2 , (14) 

It is desired to minimize the difference between U i= , and U, =2 at each point through the thickness 
in Qj. This can be accomplished by first integrating the square of the error through the thickness 
as follows. 

E 4/( eTe ) dz (15) 

a 

where a and b define the limits of integration through the thickness as indicated in the interval given 
in Eqn. 13 These integration limits must be considered carefully since the presence of surface 
bonded actuators/sensors may change the dimensions of the laminate in any of the regions. It is 
desired to find a relationship between the independent functions in Q and £1^ which minimizes the 
error in terms of the functions in the nondelaminated region to satisfy the continuity conditions. 
Therefore, derivatives of E are taken with respect to the independent functions in £> d] and are set to 
zero as follows. 


/dw . 'l 


3E 

d i — 

I dy 


(16a-f) 


Taking derivatives and rearranging the above equations leads to the following relationships which 
satisfy the continuity conditions. 
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(17a-g) 


dx 
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yi 


where 1 - 2 and the above relationships correspond to Q and n d2 , respectively. The constants otj, 
Pi, and yj are as follows. 

CXi= (- 4a 4 - 36a^ 3 b - 60a 2 b 2 - 36ab 3 - 4b 4 + 36a 2 c, 2 + 68abc, 2 + 36b 2 c, 2 + 52a 3 Cj + 228a 2 bCj + 228ab 2 c + 
52b c, - 72a-c,c, - 1 36abc,Cj - 72b 2 c,Ci - MOac^q - 140bc, 2 Ci - 156a 2 q 2 - 388abq 2 - 156b 2 q 2 + 280ac,c 2 

m?K°S C 2 C \n + < ^ ^ + , 14 ° aCi3 + 140bC|3 ' 280c ' c i 3 • 27a2d * 2 ' 51abd * 2 - 77b 2 d] 2 + 1 05aCjd , 2 + 
lOSbcjd, 2 - 1 05cj 2 d , 2 )/(3(9a 2 + 17ab + 9b 2 - 35aq - 35bq + 35c, 2 )d, 2 


(18) 


Pi - (18a 2 + 34ab + 18b 2 - 35ac, - 35bc, - 35a Cj - 35bc, + 70c,Ci)d, 2 /(2(9a 2 + 17ab + 9b 2 - 35ac, - 35bc 
35c, 2 )dj 2 ) 


Y,- (- 30a- , c, - 1 ]0a 2 bc, - 1 10ab 2 Ci - 30b 3 c, + 72a 2 c, 2 + 136abc, 2 + 72b 2 c, 2 + 30a 3 q + 1 10a 2 bc- + 1 ]Oab 2 c + 
30b-c, + 56a 2 c, Ci + 168abc |Ci + 56b 2 c, Ci - 280ac, 2 Ci - 280bc, 2 Ci - 128a 2 c, 2 - 304abc 2 - 128b 2 c 2 ‘ + 
]40aC| C) 2 + 140bC|Cj 2 + 280 C| 2 Ci 2 + 140aq 3 + 140bq 3 - 280c,q 3 - 18a 2 d, 2 - 34abd, 2 - 18b' 2 d, 2 + 70ac',d, 2 


+ 70bc,d, 2 - 70c, 2 d, 2 + 1 8a 2 dj 2 + 34abd, 2 + 18b 2 d, 2 - 35ac,d 2 
70c 1 c,d, 2 )/(2(9a 2 + 17ab + 9b 2 - 35aq - 35b Ci + 35 Ci 2 )d, 2 ) 


(19) 


35bC|d, 2 - 35ac,d, 2 - 35bc,d, 


( 20 ) 


Identical expressions are formulated which correspond to Q and <2> d3 by setting i = 3. Multiple 
delaminations can be incorporated into the developed theory by defining additional regions of 
delamination at arbitrary locations in the laminate (i.e. i = 1, 2, 3, 4, 5, 6 ...). 
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Finite element analysis 

The finite element equations are derived using the discretized form of Hamilton’s principle, 
which is stated as follows. 



where ti and t 2 are the initial and the final times, respectively and 5K C , 5U e and 5W e are the 
element variations in the kinetic, strain and potential energies, respectively. The finite element 
matrices are formulated as follows. 

>2 N e 

J £[5u c M c u e + 5u c K'u c - 5 u'(f c + F')] = 0 (22) 


where N e is the number of elements, an overdot indicates a derivative with respect to time and the 


nodal degrees of freedom for each element, u* , are specified as follows. 


U n = 


dw„ 


nT 


3x dy 


^ 4>y 


(23) 


Elemental nodes for the nondelaminated and delaminated regions are placed at the local midplanes 
of each region as shown in Fig. 6. The mass matrix, M e , is formulated as follows. 

M c = J pN' T ZN'dA' (24) 

A c 


where p is the density, N c contains the elemental shape functions and Z contains the z dependence 
of the displacement fields. The stiffness matrix, K c , includes bending and extension terms 
(subscript B) and transverse shear terms (subscript T) 

F = A B B B dA c + J B T A T BjdA c (25) 

A' A c 

where A B and A T are laminate stiffness matrices and B b and By relate the strains to the interpolated 
functions. The subscripts “B and “T” correspond to the extension/bending and transverse shear 
strains, respectively. Two force vectors are formulated for the distributed load, F' , and the 
piezoelectric forces, Fp . 

F c = jN'V(x,y)dA e (26) 

A c 
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( 27 ) 


Fp — Jn b A B ^AdA' 

A' 

The linear finite element equations of motion are obtained as follows. 

Mii + Ku = F + F P 

where the quantities M, K , F and w denote the mass and stiffness matrices, the force vector due 
to a distributed load and the nodal displacement vector, respectively. The quantity F p is the force 
vector due to the piezoelectric actuation. Bilinear shape functions are used for the inplane 
displacements and rotations (u, v, <j> x , <}) y ) while a 12 term cubic polynomial is used for the 
transverse displacements (w). The resulting four noded rectangular elements are nonconforming 
for computational efficiency and contain 28 degrees of freedom each. 

Implementation of continuity conditions 

The continuity conditions between the nondelaminated and delaminated regions need to be 
included in the finite element model by constraining the nodal degrees of freedom at the interface of 
the nondelaminated and the delaminated regions. Applying these constraints directly to Eqn. 28 
leads to a nonsymmetric system of equations which results in undesirable consequences in the 
solution technique. A different approach must be used. The continuity conditions presented in 
Eqns. 17a-g for the regions and Q d2 , as well as similar conditions for the regions Q and Q d3 , 

are applied to the finite element degrees of freedom at the interface of the nondelaminated and 
delaminated regions (E>i) presented in matrix form as follows. 


Ru = 0 


(29) 


where u, which is a subset of the u, contains the degrees of freedom in £2] and is partitioned as 


follows 


U 


u i=l 
u i-2 
u i = 3 


(30) 


and u, contains the nodal displacements in the i-th region (i.e. i = 1 refers to Q as before) which 
are also contained in as follows. A similar relationship is defined for u. The matrix R has the 
following form. 
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R = 


R: 


i = 2 "I O' 

.R, = 3 0 -I 


(31) 


where the first column contains the submatrices Rj =3 and Rj =3 which in turn contain factors in 
terms of a i( bj and h. The exact formulation for these submatrices will be presented in the final 
paper. The Lagrangian for the discretized system is formulated as follows 

L = ^(u t Mu - u t Ku + u t F + u T F p ) + 1 p(u T R T RG + u T R T Ru) ^ 2 ) 


where the first four terms are the kinetic, strain and potential energies due to an applied load and 
piezoelectric actuation, respectively. The last two terms are penalty terms for the displacements 
and velocities multiplied by a penalty factor, p, to satisfy the continuity conditions at £2j. 
Differentiation of L with respect to the nodal degrees of freedom using the Lagrange technique 
allows the equations of motion to be reformulated to include the constraint terms while retaining 
symmetry. The resulting augmented equations are as follows. 


[M + pP]ii + [K + pP]u = F + F P 
where P is the penalty matrix defined as follows. 

P = R T R 

which is expanded to correspond to the global degrees of freedom (u). 


(33) 


(34) 


Preliminary Results 

Static response 

Some numerical and experimental results are presented for a cantilever Gr/Ep composite plate 
with stacking sequence [0°/9070°] 6 and dimensions L = 25cm, W = 10cm, and h = 2.286mm. The 
thickness of each ply is 0.127mm which is typically used by industry. Plies numbered 13-16 are 
replaced with a single PZT piezoelectric layer of thickness 0.381mm which represents a 
commercially available piezoelectric material thickness. Material properties are presented in Table 
1. The plate mesh has 25x10 elements and 2002 degress of freedom. The undeformed mesh with 
no delaminations is presented in Fig. 7. The piezoelectric layer is energized with a constant voltage 
of 150 volts. Both longitudinal bending and camber deflections are present due the actuation as 
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shown in Fig. 8. This is because the piezoelectric layer is offset from the neutral axis of the plate 
which produces bending moments about the X and Y axes. The finite element results correlate 
well with sample static experimental tests as shown in Fig. 9. In this figure, the transverse 
displacement is denoted (w*) and is normalized with respect to its maximum value. The 
lengthwise coordinate is denoted (x*) and is normalized to the length of the plate . 

Next, the effects due to delamination are investigated using the above smart plate properties. 
Delamination is introduced at the interface between the composite substructure and the piezoelectric 
layer (plies 12 and 13) at the free end of the plate through the entire width ( 20 cm < x < 25 cm, 0 
cm < y < 10 cm). The mesh for the composite substructure, located at the midplane of the 
nondelaminated portion of the plate, comprises 20x10 elements. Two additional 5x10 meshes are 
generated above and below the delamination at the midplanes of each respective region. These 
meshes are offset from the midplane of the undelaminated portion of the plate as shown in Fig. 10. 
It is important to note that Figs. 10 and 1 1 are not to scale so that the different meshes may be 
easily identified. Continuity at the interface of each of these delaminated regions and the composite 
substructure region is obtained by implementing the continuity conditions and penalty approach 
discussed earlier. The total number of degrees of freedom for the delaminated plate is 2541. 
Again, the piezoelectric layer is energized with a constant voltage of 150 volts. The static 
deflection in this case (Fig. 1 1) is significantly different from the first case (Fig. 8) where actuation 
produced bending which resulted in negative transverse displacements. The nondelaminated 
portion of the plate bends downward as before, as does the region below the delamination. 
However, the region above the delamination bends in a positive transverse direction. This is 
because the offset of the piezoelectric layer from the local midplane in the delaminated region is 
opposite the offset in the nondelaminated region. This results in a reversal of the bending moment 
due to the piezoelectric actuation and a positive transverse displacement. These preliminary results 
indicate that a thorough understanding of the effects due to delamination is critical to designing 
structures with smart composite materials. The final paper will present additional results to 



investigate the effect of delamination on the static and dynamic response of a smart composite 
laminate. 

Dynamic response 

In this section, results are presented to demonstrate how the developed theory can be used to 
investigate issues related to active control and delamination. The test article is a composite 
cantilever plate with dimensions W = 5.08 cm, L = 45.72 cm, and thickness 1.16 mm. The 
material properties are again those of Gr/Ep (Table 1) and the stacking sequence is [0°/90°]4s. Two 
piezoelectric actuators with dimensions 2.54 cm (width) x 10.16 cm (length) x 0.762 mm 
(thickness) are bonded side by side on top of the laminate near the root of the plate as shown in 
Fig. 12. The actuators are linked electrically to function as a single unit. The piezoelectric sensors 
have the same width and height as the actuators. However, the thickness of the sensors is 0.381 
mm. They are placed opposite the actuators on the bottom of the laminate and also function as a 
single unit. The first natural frequency of the laminate, which corresponds to the first bending 
mode, is 0)] = 9.1887 Hz. A damping ratio of = 0.02 is assumed to represent the passive 
damping capability already present in the composite material. The laminate is initially loaded with a 
static distributed force of q = 30 kg/m 2 which is suddenly released at time t = 0. The sensor and 
actuator are coupled to a velocity feedback control loop which is designed to reduce vibratory 
amplitudes during the transient response as quickly as possible. Practical limitations dictate that the 
voltage supplied to the actuator must not exceed 100 volts. This corresponds to an electric field of 
138 V/mm which is well below the maximum allowable electric field of the piezoelectric material. 
Five modes are included in the dynamic response analysis. 

The test article is first examined with the control system turned off. The tip of the laminate is 
displaced 1 mm due the static load and is then released. The decay envelope due to passive 
damping alone is presented in Fig. 13. In this figure, the transverse displacement at the tip of the 
laminate (w) is presented for both the passive case and the active case which will be discussed 
shortly. The 1 % settling time (t s ), which is the time when the vibrational amplitude become \ % of 
the initial amplitude, is determined from the following equation. 
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( 34 ) 


_ logc(OOl) 

1 r. ~ 

^i“i 

The ]% settling time for the passive damping case is found to be 3.732 seconds. 

Next, the velocity feedback control system is turned on and the gain is selected to appropriately 
energize the actuators without exceeding the maximum allowable voltage. The voltage supplied to 
the actuators during the transient response is shown in Fig. 14 and does not exceed 100 volts. The 
active controls significantly increase the damping ratio of the composite laminate from 0.02 to 
0.0513. This reduces the settling time by 61% to 1,455 seconds as shown in Fig. 15, The 
transient response for both passive and active cases are presented in Fig. 13 Thus, the 
piezoelectric sensor/ac.uaior pair has a significant effect on increasing the damping of the laminate. 

Now the effect of delamination on the control system is investigated. Deiamination is 
introduced in the above test article at the tip of the actuator nearest the fixed end and extends 
throughout the width of the actuator at the interface of the bonded actuator and the composite 
substructure as shown in Fig. 16 The length of the delamination (L d ) is indicated by the 
nondimensional parameter P where P = L„/L and L is the length of the laminate. As the length of 
the delamination grows, the control authority is reduced since the actuator is less able to transfer 
the stresses, which are induced from the piezoelectric actuation, to the substructure. This reduction 
in control authority is indicated in Fig. 17 where the settling time is presented for a range of 
delamination lengths. The current model allows for both slipping and separation of delaminated 
layer (in this case the bonded actuator) from the composite substructure. Separation of the actuator 
a. p = 1 1 % is observed the sixth mode, which is the fourth bending mode, a, a frequency of e) 6 
= 209.8 Hz as shown in Fig 1 8. Separation of the delaminated actuator layer ,n one of the lower 

modes is unexpected and may significantly inDuence a more complex dynamic response, such as 
the case of a forced response. 

Future work 

Research plans that are either already under progress or are soon to be implemented are as 
follows. 





• Complete mathematical investigation to study the use of smart materials for vibration control. 

• Determine the effect of delamination on the authority of the control system. 

• Develop and implement experimental procedure to validate proposed mathematical model. 
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Table 1. Material properties. 



Gr/Ep 

PZT 

Ei (GPa) 

143 

63.0 

E 2 (GPa) 

9.70 

63.0 

V 1 2 

0.300 

0.300 

G 1 2 > Gj 3 (GPa) 

6.00 

24.2 

G 23 (GPa) 

2.50 

24.2 

p (xlO 3 Kg/m 2 ) 

1.39 

7.60 

d 3 ] ( xlCH 2 m/V) 

- 

253 
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Figure 1 Composite laminate cross section 
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Figure 2 Smart composite plate incorporating piezoelectric layers. 
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Figure 3 Laminate cross section. 
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Figure 4 Local region plate geometry. 




Figure 6 Finite element discretization. 
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Figure 10 Undeformed mesh including delamination. 
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Figure 1 1 Static deflection due to piezo actuation including delamination. 
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Figure 12 Gr/Ep laminate with bonded piezo sensors/actuators 
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Figure 1 3 Passive / active transient response. 




Figure 15 Decrease in settling time due to active control. 
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Figure 16 Composite laminate with delaminated actuators. 



Figure 17 Increase in settling time due to delamination. 
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Figure 18 Separation of delaminated layer in mode 6 
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